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1. Introduction 

No known experimental constraint firmly excludes the possibility that Kaluza Klein 
(KK) excitations of the graviton propagating in 5 > 2 large extra dimensions will 
affect future particle physics experiments |1|, 0^. After removing all non propagating 
degrees of freedom by a suitable choice of coordinates, many authors computed the 
signals of KK graviton emission at tree level 0, ^1, § • Some authors also considered 
1-loop effects ^, H, ^ [l^: since they affect observables measured with higher 
precision, they can compete with tree level effects. The result was not the expected 
one. Consider for example the graviton correction to the Higgs mass. At first sight 
one would estimate it as 

where M^, M4 are the gravitational scales of the D-dimensional and 4-dimensional 
theory respectively, and A ~ Md parameterizes the unknown quantum gravity ul- 
traviolet cutoff. At a closer look @, 0, ^, |^ the effect seems to be much larger. To 
understand that, consider the propagator for the physical J = 2 nth KK graviton 
with mass rrin and 4-momentum |^ 

where t^^, = rj^^, — k^ky/m\. If the terms enhanced by powers of /c/m„ were to fully 
contribute to quantum corrections, the k factors would give a highly ultraviolet (UV) 
divergent loop effect. More importantly, when 5 < 4 the l/m„ factors would also 
give a strong infrared (IR) enhancement of the sum over KK modes. At the end 
the correction would be a factor (M^ji?)^"^ larger than the naive one in eq. ( |1 . 1| ) , 
where R is the size of the extra dimensions. This kind of behavior, indeed observed 
in 0, H, 1^, would exclude the possibility that 5 < 4 large extra dimensions (i.e. 
R ^ I/Md) solve the hierarchy problem. 

The above argument on the fate of the k/mn terms must however be wrong, 
and for a very simple reason. Indeed one could choose to fix the D-dimensional 
reparametrization invariance by the de Bonder gauge choice, in which the graviton 
propagator contains no k/mn terms. This is in complete analogy with the case of 
a massive vector boson, where the propagator contains k/mn terms in the unitary 
gauge, while no such term is present in the Feynman gauge. Therefore fc/m„ terms 
cannot affect gauge-invariant physical observables. This suggests that there must be 
something missing or incorrect in the computations so far performed. 

^The case 5 = 2 is excluded by bounds on emission of KK modes with a small mass ^ 100 MeV 
in supernovae ft], ||, if the extra dimensions are flat. In principle, one could save collider signals 
(due to heavy KK modes) by assuming that the compact dimensions are curved on length scales 
> (100MeV)"\ so that the light KK are hfted. 
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The purpose of the present paper is to devise all the elements that are needed 
for a fully consistent computation. The guideline is to respect the full D-dimensional 
general coordinate covariance. First of all it is crucial to fix the gauge by the Faddeev- 
Popov procedure and to choose a covariant regulator. If the regulation of the loop 
integral is not performed with the due care, spurious UV and IR divergences can 
appear. Secondly, one has to remember that the position of the brane depends on the 
system of coordinates, and therefore brane fluctuations (branons) must be taken into 
account in order to respect general covariance. Finally one must carefully identify 
which are the physical observables in the presence of gravity: misidentification of the 
true observables can yield spurious gauge dependence and IR divergences. 

One of our results is that all the puzzling effects found in the existing literature 
cancel out in a fully consistent calculation when one computes physical observables. 
For example the Higgs mass term and the oblique S, T, U parameters ||rT| are not 



physical observables (except in particular cases). So they receive gauge dependent 
quantum gravity corrections, which in some cases are even enhanced by powers of 
RMd. These infrared pathologies, which would invalidate perturbation theory (for 
instance RMd ~ 10^^ if i? ~ mm and ~ TeV), are absent in the corrections that 
affect the corresponding physical observables, the pole higgs mass and the €1,62,63 
parameters ]12| . 



In our study we treat quantum gravity and the brane by the method of effective 
field theory (EFT) [|13|, [l^- We do so in the absence of a realistic fundamental de- 



scription of the SM on a brane^. The effective Lagrangian summarizes all our low 
energy knowledge of gravitational interactions with SM particles. By our method we 
could perform a fully consistent computation of the 1-loop quantum gravity correc- 
tions to electroweak precision observables. However the dominant effects are strictly 
speaking uncalculable, as they are saturated in the UV where we loose control of the 
theory. We can only parameterize these effects in terms of a UV cutoff A^. The calcu- 
lable piece is the one saturated in the infrared, but this is only of order {Mz/MdY^^- 
Therefore, introducing a UV cut-off A, we will only compute a particular combina- 
tion of observables, which is affected by just a few simple Feynman diagrams. For 
the full set of observables we will limit ourselves to a qualitative discussion. 

While the discussion of the phenomenology is somewhat limited by the powerlike 
UV divergences, we stress that the main goal of the present paper is conceptual. In 
this respect the most important (and new) result is that brane motions have to be 



^Interesting attempts based on D-brane intersections JTsI give 'semi-realistic' models with extra 
charged matter with respect to the SM. The stability of these configurations is an open question. 

■^A string model could provide a physical realization of this cut off. However at the level of the 
present model building technology there are many free parameters specifying the moduli and the 
brane configuration ||l^, [l^, |l^. Therefore, even if we were able to reproduce the SM, the predictive 
power on quantum corrections (for example on the muon g — 2) would probably be limited. Of 
course it would still be important to have one such model. Indeed it would also be interesting to 
have a field theoretic brane model in the spirit of HS] . 
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properly taken into account. In order to understand this issue better we have also 
considered the case of a brane living at an orbifold fixed point, for which the branons 
are projected out. In this case gauge independence of observables is met through 
tadpole diagrams specific of orbifold compactifications, rather than by branon loops. 
The technology developed in this paper may prove useful in future work. One possible 
application is the brane to brane mediation of supersymmetry breaking through bulk 
gravity at 1-loop. This effect is computable and represents the leading correction to 
anomaly mediated soft terms: depending on its sign it may cure the tachyon problem 
of anomaly mediation. 

The paper is structured as follows. In section 2 we discuss our Lagrangian 
and the effective field theory philosophy. We also introduce various gauge fixing 
conditions for the gravitational field and explain the role of the branons. In section 
3 we calculate the corrections to the masses of scalars and vectors, on and off the 
brane. We explain what gauge independence means in quantum gravity, and show 
that physical quantities are gauge independent. We also give the example of a brane 
at an orbifold fixed point, for which branons are not needed. In section 4 we derive 
experimental bounds on low energy quantum gravity from precision measurements 
and from the anomalous magnetic moment of the /i. In section 5 we summarize. 
Finally in the appendices we describe how to derive graviton-matter vertices and 
collect our results for the corrections to brane observables. 

2. Effective Lagrangians for gravitons and branes 
2.1 Pure gravity 

We study gravity in x where the extra dimension is a compact manifold of 
dimension S. Not knowing which manifold is of physical interest (if any), we consider 
the simplest one: a 5-torus with a single radius R and volume V — {27rRy. We 
perturbatively expand the classical Einstein-Hilbert action around the flat metric 
Qmn — Vmn + i^hMN: Vmn — (+1, —1, —1, —1, ■ ■ ■) in terms of the graviton field 

hMN' 

]lfD-2 p 

= ^ y d^'X ( - h^^'uhMN + hnh - 2h^''dMdNh + 2h^'^dRdsh^M) + 

(2.1) 

where D — d + 5, h = and we used tjmn to raise and lower indices. We use upper 
(lower) case latin letters for D-dimensional (extra-dimensional) indices and Greek 
letters for d-dimensional indices; in particular we decompose the £)-dimensional co- 
ordinates as — {x^, y^). We do not fix d = 4 since we will use dimensional 
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regularization. Following the notation of ref. [^] we have defined 

/t^ = 4M|,-^ , M^^ = VM^-^ = R^Mg-^ (2.2) 

Md is the effective reduced Planck mass as measured by a d-dimensional observer, 
Mr) is the corresponding parameter in D dimension and Md is defined by ( |2.2| ). 
With this convention the equations of motion read: 



1 „ 1 



Rmn — -^dMNR — — -,g_2 Tmn (2.3) 



Before inverting the quadratic term in eq. ( p.l| ) to obtain the propagators, one must 
fix the reparametrization invariance; we follow the Faddeev-Popov procedure and 
introduce a set of ^-gauges by adding to the Lagrangian the gauge-fixing term £gf = 
—F'^/^, where 

Fn = d>'{h,^ - ^v,Nh) + ^d\KN - Y^mNh). (2.4) 

This particular choice breaks the D-dimensional Lorentz symmetry of the fiat back- 
ground metric for generic values of ^ and interpolates between the usual de Bonder 
and unitary gauge, obtained respectively in the limit ^ —>■ 1, oo. The functional inte- 
gral gets multiplied by the Faddeev-Popov determinant, exponentiated in the usual 
way by introducing 'ghost' fields rjM, Vm- 

£ghost = / d^X d^X' miX)^^^ VMiX') (2.5) 

where A is the gauge parameter for reparametrizations. 

The kinetic term for the graviton field is a (messy) 3x3 matrix which mixes 
tensor /i^j,, vector /i^j and scalar hij modes. Since interactions are more easily written 
in terms of the /i^j^, /i^j and hij components of the Z)-dimensional graviton field Hmn, 
it is more convenient to write the propagator in this basis rather than in the gauge- 
dependent mass eigenstate basis. For example matter fields confined on a straight 
(i-dimensional brane at leading order couple only to the tensorial hf^i, mode. 

By decomposing the graviton field Hmn in its Fourier harmonics 

hMNix, 2/) = i E ^MNi^^ e'""'^^ (2.6) 

and integrating the Einstein-Hilbert Lagrangian over the extra-coordinates, one ob- 
tains the (i-dimensional Lagrangian for KK modes. 

Notice that d'^hn^] can be interpreted as Goldstone bosons eaten in a gravitational 
Higgs mechanism to form massive tensors and vectors. We are classifying particles 
by the d-dimensional Poincare group. By this interpretation, eq. (p.4|) is the analogue 
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of 't Hooft's ^ gauge in spontaneously broken gauge theories. For ^ oo we get the 
unitary gauge 0, |[ in which only the physical degrees of freedom propagate. In this 
limit the ^ gauge propagator for the modes Hmn = {h^^, hij, h^j) simplifies to 



^ ^ • 1 l6r, 



An) An') _ lOn-n' 
'''MN'''M'N' ~ ~' 



2 (P - ml) 

2 pf P P A_ p p _ 2 p p fl 

-Pn'V 



(2.7) 



where /c is the rf-dimensional momentum 



mi 



and = / is the mass squared for the nth KK excitation, having defined 
71^ = —n.injT]'^^ = niTijdij. In appendix A we derive this propagator by working in 
the unitary gauge with physical fields. As usual 'Goldstone' bosons and 'ghosts' get 
infinitely massive when ^ — > oo but they do not decouple: loop corrections computed 
in the unitary gauge = oo) by propagating only the physical fields are different 
from the limit ,^ — cxd of loop effects computed in a ^ gauge ||20|. Of course the 
mismatch disappears in physical quantities. 

For ^ = 1 we get instead the de Bonder gauge, where the propagator has the 
covariant form: 



I 1 i 2 

hMNhM'N' = „ J- {vmm'Vnn' + Vmn'Vnm' — 7^; ^VmnVm'N') (2.9) 

ZK ^ U — I 

where K is the D-dimensional momentum. In matrix notation, for the single KK 
mode: 



7 (n) 7 (n') 1 i^n,—n' 

'^MN'^M'N' ~ ~ 



2 {k"^ — ml 

2 X ^ X X I X X . 2 



-5iiirjf,^> 

(2.10) 

We have thus shown that the propagator in the de Bonder and unitary gauges has 
the same form up to longitudinal k/mn terms. For compactness, we do not write 
explicitly the propagator in a generic ^-gauge. 
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2.2 Gravity and branes 



We want to study quantum gravity corrections to the physical observables of a field 
theory living on a (i-dimensional brane in a (5-dimensional compact space Ai . At the 
end we will identify the brane theory with the Standard Model. The gravitational 
Lagrangian has been discussed in the previous subsection. We now discuss the brane 
Lagrangian. 

We use an effective field theory (EFT) approach where the fundamental descrip- 
tion of the particles and of the brane is not specified [|1^]. In the regime of validity 



of EFT, the particles are treated as point-like and the brane is treated as infinitely 
thin in the extra dimensions. This requires a little explanation. If p is the brane 
true transverse size, our EFT is only valid at energy scales ^ 1/p. The brane also 
generally has a finite tension t = f^. This gives rise to a gravitational field behaving 
like f'^/M'^^-'^r^-'^ = {rc/rY''^ at a distance r far away from the brane in the extra 
space. We focus on 5 > 2 (for 6 = 2 the background is locally fiat with a conical 
singularity at the brane position). The gravitational radius controls the distance 
at which the geometry is curved. One can then think of different possibilities for 
the brane structure. If p > tq the brane is similar to a big star where the geometry 
nowhere strongly deviates from the flat, and 1/p truly represents the UV cut-off of 
our EFT. On the other hand for p ^ it is the gravitational radius that sets the 
UV cut-off. Physics at energies > l/r^ would probe the gravitational structure of the 
brane, which is non-universal and model dependent. One example is a black brane 
where at r ~ a black- hole horizon is present. Another different example is given 



by the solution studied in [£T|, where there is no horizon and a naked singularity 
is avoided by a flnite brane size. In the latter case the coupling of bulk gravitons 
to the brane is dramatically changed at energies > l/r^. As we are only interested 
in universal features we will assume that the UV cut off Auv that limits the use of 
our EFT is bounded by min(l/p, l/rc). In the regime of validity of EFT we can 
treat the background metric as approximately flat and treat the effects of the brane 
tension as perturbations. Notice indeed that at energy E the latter are controlled 
by the small parameter E^-"^ f/M^^^~^ = {EtgY^'^ < (Auyrc)^"^ < 1 . 

Two possibilities are given: either the brane can freely move in the bulk or sit 
at a flxed point, if the compact space M. has any. Let us consider the former case 
flrst. The immersion of the brane in the D-dimensional space is parameterized by 
D functions X^{z), where are the d local coordinates on the brane. The brane 
action must be invariant under both D- dimensional coordinate changes (under which 
X*^ transform and are unchanged) and under reparametrizations of the brane 
coordinates z^. An invariant brane action can be built using the induced metric 

= ' 9mn{X{z)). (2.11) 
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Since g""^!^ is a scalar under D- dimensional reparametrizations, we only need to re- 
spect brane reparametrizations by the use of g^^^ itself. The description of the brane 
position by the X*^(z) is of course redundant. We can eliminate this redundancy 
by using the remaining gauge freedom represented by brane reparametrizations. We 
stress that we cannot use D-dimensional diffeomorphisms for which the gauge has 
been completely fixed in the previous section. A convenient choice of brane coordi- 
nates is = z^, ?/* = ^^{x^). This choice completely fixes brane reparametrizations 
without the need of introducing additional ghost fields (the ghost determinant is 
trivial) [jl4| . We call the i^* hranons. 

As we said the branons cannot be thrown away because we have already com- 
pletely fixed the D-dimensional reparametrization gauge invariance. However in the 
previous section one could have chosen a different class of coordinate gauges, one 
in which the brane always sits at a given point in M.. This different choice would 
explicitly break translation invariance in the extra dimensions. What becomes of the 
branons in these different gauges? They are still there but as longitudinal modes of a 
combination of graviphotons: the branons can indeed be interpreted as the Goldstone 



bosons of broken translation invariance in the extra dimensions |]T4| . We find it more 
convenient to gauge fix the graviton in the more standard way and keep the branons. 
Notice that, consistently with their Goldstone character, in the limit in which grav- 
ity decouples {Mo 06) the branons survive. Their physical effects can therefore 
be studied independently of gravity |^. Quantum fluctuations of the branons are 
controlled by 1/r (the analogue of 1//^ for pions) and become non-perturbative at 
an energy E > a/Itt/ {E > Anfj^ for pions). Therefore the tension r sets another 
sure upper bound on the regime of applicability of EFT. 
In terms of the branons ^* the induced metric is 

gfu = - 9,^9uj9'' + {D,^i) {D^Qg'' = + V (2-12) 

where D^^i = d^^i+g^i and the metric guN is evaluated at the brane location = 
For 1-loop computations we need h^^ up to quadratic order in ^ 

V = '^V + {d^^i){duC) + <Cd^K, + h^^d^C + h,d^C) + ■■■ (2.13) 

where now h is the graviton field evaluated at the brane rest position = 0. The 
brane Lagrangian is given by 



-TJdetg-^^^ + CsM + 



(2.14) 



where £sm is the covariant brane Lagrangian (that we will identify with the SM 
Lagrangian), while the dots indicate all terms involving higher derivatives, the Rie- 



mann tensor for the induced metric |]23| or the extrinsic curvature. By expanding 
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the tension term up to quadratic order in the branons and gravitons we get 



r 
2 

(2.15) 



which shows a mixing between ^ and h (see appendix B for the definition of the tensor 
B). As we will discuss shortly, in order to consistently compute virtual graviton 
effects this mixing has to be taken into account. Notice also that there is a linear 
term in h, since the massive brane is a source of gravity. We will comment below 
about when and how can this term be neglected. The interaction of gravitons and 
branons with SM fields is encoded in the covariant dependence of £sm on the induced 
metric. At quadratic order we have 

C^M = LsM + nL^%, + K^L^^'^'^'hp.h^,, + ■■■ {T^, = -2^) (2.16) 

where the explicit formulae are given in appendix 0. 

2.3 Gravity and vector bosons 



In section ^3] we have described the gauge fixing procedure for a theory of pure grav- 
ity. If gauge fields Am are present, the Lagrangian has both internal and gravitational 
gauge invariance, which can be fixed through a delta functional 6{F{h, A)) 

F{h,A) = [Mh,A),Mh,A)] 

in the functional integral imposing fi{A,h) = 0, f2{A,h) = 0. This is equivalent 
to adding the gauge fixing term £gf = —fiH — /l/^C i^i the Lagrangian and the 
Faddeev- Popov determinant in the functional integral 

where Ai, A2 are the gauge parameters for diffeomorphisms and internal gauge trans- 
formations respectively. 6/2/6X1 is generically non zero because vector bosons are 
'charged' under gravity. However the graviton field is neutral under charge transfor- 
mations, so that for a reasonable gauge-fixing function /i which doesn't involve the 



vector bosons (in particular for /i as in eq. ( |2.4]) ) the determinant factorizes: 



det ^ = det |A . det 1^ (2.18) 
OA oAi 0A2 

and the two factors can be exponentiated separately in the usual way. Notice that 
it is convenient to choose a non covariant gauge-fixing /2 for the photons in order to 
avoid additional couplings with gravitons. (A non-covariant /2 should not cause any 
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panic: reparametrizations are already broken by the gravitational gauge fixing fi). 
We have explicitly checked that the simple gauge fixing 

^GF = -^(9mAV^)' (2.19) 

gives the same results as other more involved choices. 

In a theory with vectors that acquire mass My through the Higgs mechanism, as 
in the Standard Model, the gauge-fixing term will contain the Goldstone bosons (pc 
field as well. Even with the simple gauge fixing 

Cgf = -^{OmA^t]^''' - MyC<pG? (2.20) 

there is a cubic vector-Goldstone-graviton interaction: in a generic metric the gauge 
fixing does not fully cancel the kinetic mixing between the Goldstones and the vector 
{MyAMdN<t>G9^''^ ^/g) present in the Lagrangian. Such gauge fixing can be easily 
adapted to vector bosons confined on a brane. 

2.4 Gravity and fermions 

Finally, we sketch how to extend our analysis to the important case of fermions. It is 
well known that GL(D) does not admit spinor representations and in order to deal 
with fermions, we need some extra structure: the vierbein E\j and its inverse E^^ 
defined by 

Qmn = Vab Ej^jE^j^ Ej^fE^^ = 6^ . (2.21) 

Where capital letters from the beginning of the latin alphabet A, B,C, . . . denote D- 
dimensional Lorentz indices. The vierbein basis definition introduces an additional 
gauge symmetry, besides diffeomorphisms, due to the freedom in ( |2.21| ) to rotate E 
acting with a local SO{D — 1, 1) transformation. In absence of torsion, the compati- 
bility condition between the metric and the connection u, allows to express the latter 
in terms of the vierbein E. Then, once the vierbein is defined, the introduction of 
spinors is rather straightforward (see for instance [0]), a collection of the relevant 
formulae can be found in appendix B. Around a fiat background we can parametrize 
the vierbein as E^.j- = + kB\j. In terms of the quantum field B\j the metric 
fiuctuation is then 

h,,N = B,,, + 5,,, + kBo^^B^^. (2.22) 

where Bmn = VmaB^ and similarly all indices are raised and lowered by the 
Minkowski metric tjab- The gravitational action, when expressed in terms of E 
(or equivalently in terms of B), is invariant under the infinitesimal local Lorentz 
transformation 

6Bij = K-'ni{x) {6fj + kbIj) , nAB{x) = -^BAix). (2.23) 
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A convenient gauge choice is p5| 



^M^. - = 0, (2.24) 



The great advantage of (|2.24| ) is that Lorentz ghosts are absent |2^ and that it makes 



possible the ehmination of the vierbein fields, order by order in k, in favor of the 
quantum metric h |2^. Indeed, in the gauge ( |2.24| ) one can easily express B in terms 
of h by solving ( [2.22| ) 

= Ik.. - h\^h,, + O(k^). (2.25) 

As a result, even when fermions are present, at the perturbative level, the quantum 
fluctuations of the geometry are encoded in h and our formalism can be applied 
without modifications. 

A similar procedure can be applied to fermions living on a (c? — l)-brane. These 
are spinors of SO{d — 1, 1) and in order to write an invariant Lagrangian one needs 
the induced vierbein on the brane, which is now a d x d matrix e";,. In what follows 
we indicate by lower-case latin letters a,b,c, . . . the (i-dimensional Lorentz indices. 
In ref. |jT^ it was shown how to construct e"^ out of and of the brane immersion 



X (z). Basically it has the form 

e;(^) = (Xiz)) d,X'\z) (2.26) 

where is a SO(D — 1, 1) rotation matrix which depends on E\j and X^ [z). Under 
a SO(D — 1, 1) rotation E\j — > VL{X)^E\j, the induced vierbein undergoes a S0(o? — 
1, 1) rotation e"^ — > uj{z)le^^. By fixing the brane reprametrizations keeping just 
the branons (as done in the previous section) and by fixing D-dimensional Lorentz 
transformations as shown in this section, is written as a function of ^^{x) and 
Hmn- However it is a fairly complicated expression. Calculations can be simplified 
by using the local Lorentz symmetry — >• ujl{x)e^^ to rotate the induced vierbein to 
a more convenient form (fermions rotate ip"" — >■ u'^ip^ by the spinorial representation 
u'p). Precisely as we did with E"\.]- it is useful to rotate e"^ to a symmetric matrix 

= 5; + 6; v.ab\ = v,abl (2.27) 

from which by using (^"jf = e" e"'' and eq. (|2.12|) we obtain the analogue of eq. ( p.25| ) 



b^u = - + 0{V). (2.28) 

Z o 

3. Loop corrections to brane observables 

We now have all the ingredients to perform some illustrative computations. We will 
focus on the one loop correction to the masses of scalars and vectors living on the 
brane. 
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In order to do a meaningful computation we must employ a regularization that 
respects D- dimensional reparametrization invariance. The result will depend on the 
choice of the regulator. The simplest thing could be cutting the loop integrals at 
A. Since this is not an invariant regulator, we would get a meaningless result that 
also depends on the choice of the loop integration momentum. A better possibility 
consists in dividing all graviton propagators by some power of (1 — p^/A^). It is 
possible to obtain this Pauli-Villars (PV) regulator in a covariant way by adding to 
the action suitable higher derivative reparametrization invariant terms involving just 
the metric. We will instead employ the standard extension of dimensional regular- 
ization to the case in which both continuum and discrete momentum are involved 
(see the appendices for details). Of course by this method we are only sensitive to 
the "physical" logarithmic divergences, while all power divergences are automatically 
removed. Nonetheless from our results it will be clear that by choosing a regulator 
sensitive to power divergences (like PV) for the sum over KK we would still not have 
the AR terms of ref.s |^, H, in physical quantities. 

3.1 Brane in a torus 

Consider now the one-loop graviton correction to the pole mass mo of a minimally 
coupled scalar living on a straight brane located at the point = of a torus 
T^. As we explained in section |2.2| , in the regime of validity of EFT {E < I/tc) 
the brane tension can be treated as a perturbation in the gravitational dynamics. 
Therefore it makes sense to expand the corrections to our observables in a power 
series in r. Let us focus on the lowest order effects, i.e. those that go like r° ^. The 
diagrams that contribute at order r° are shown in fig. 1. Notice that diagrams (d) 
and (e) also involve branons: in these diagrams the from branon propagation 
is compensated by the in the graviton-branon mixing insertion. Notice also that 
the tadpole diagram (c) gives no contribution. Due to momentum conservation in 
the extra dimensions (valid at zeroth order in r) only the zero modes mediate this 
tadpole, these are the 4d graviton and the radion. Whatever mechanism stabilizes the 
radion giving also a vanishing effective 4d cosmological constant generates a tadpole 
that cancels (c) exactly at the minimum of the radion potential. Of course exact 
cancellation of the 4d cosmological constant requires the usual fine tuning. Now, the 
genuine graviton diagrams (a) and (b) give a correction 

Sml{a + h) = ^J2 f -2L^^L'''^ + 4^L^'^^'^|0(p))(/^jr/i(;")) (3.1) 

■^Notice that there are also corrections from pure branon exchange, which go hke inverse powers 
of r and which persist when gravity is turned off. The lowest, physically meaningful correction of 
this type to the scalar mass comes at two loops and goes like Smo/nio ~ ttiq/t^. The 1/r effects 
can be bigger than the gravitational ones we study, but they are physically independent . Thus 
it makes sense to focus only on the latter. 
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Figure 1: One-loop gravitational corrections to the pole mass of a scalar on a brane 
from gravitons (diagrams a,h,c) and from graviton/hranon mixing (diagrams d,e) at zeroth 
order in the brane tension. The mass of a vector particle on a brane also gets corrections 
from vector- Goldstone-graviton vertices (diagram a'). Gravitons (branons, scalars, vectors, 
Goldstones) are drawn as pig-tail (dot-dashed, dashed, wavy, dashed-wavy) lines. 

where = ra^ is the squared momentum of the on-shell scalars (mo is the tree level 
mass). The branon contributions (d) and (e) give 

+ = '{^(-p)\T,,\<p{p))Y, [ K{k) (3.2) 

where (0(— p)|T^^|0(p)) = 2ml Fn{k) represents the contribution of the loop 
in fig. (ld,e). The sum of all contributions in the interpolating ^ gauge defined in 
section 2 gives 

K = ^2M',-\^{d + 6-2) ^ {^^^' ^^"^'^"^'^ - - 2)(5 - 4)Ao(mg) + 
4d [{d - 2)ml -2{d + 6- 3)ml] Bo{ml, ml, ml) + d6{d - 2)mlBi{ml, ml, ml) } 

(3.3) 

/(e, d, 6) =m - ^f'' [2 - 6{2^ - l)(e - 1) - e(10 - 3^; + S^d - 4))] + 
2^''/' [3d6 + d^{6 - 2) + d\^ + 1) + 4 - 12^ + 26{^ + 6-4)] + 
2{d - 2)[ - 6 + 2d{d + 6 - 2)] 

where the Passarino-Veltman functions Aq, Bq, Bi are defined in appendix D, where 
we describe how the cutoff-independent contribution can be extracted. The mass 
correction is multiplicative as expected for a minimally coupled scalar: for vanish- 
ing tree level mass the scalar is derivatively coupled to gravity. Although in this 
expression all the terms enhanced by l/m^ found in P, H, ^ cancel out mode by 
mode, we do not obtain a gauge-independent result. However the ^ dependent term 
in 6ml/ ml, only depends on and R (and the UV cut-off Ayy if dimensional 
regularization is abandoned) but not on ml itself. So it looks like a universal effect. 
Indeed one finds the same gauge dependent piece in the correction to the mass of a 
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vector particle. In short the correction to the pole mass of a spin s = {0, 1} particle 
on the brane can be written as 



6m'i = 2mlG{MDR) + m; 



(3.4) 



where G is the only gauge dependent factor, while Ao,Ai are gauge-independent 
(/2 is the renormalization scale). Explicit expressions for these functions can be 
found in appendix |C[ Similarly we find that for a localized photon the gravitational 
correction to the electric charge e has a gauge-dependent factor equal to [e]G, where 
[e] = 2 — d/2 denotes the dimension of the electric charge in d dimensions. The 
moral of these results is that the gauge dependence can be reabsorbed by changing 
the normalization of the graviton field Qmn- In niore physical terms, gauge dependent 
terms amount to a change of the mass unit: all the dimensionless quantities that we 
have computed (like mo/mi) are gauge independent. In the presence of gravity only 
dimensionless quantities are real observables, as they are invariant under rescaling 
of the metric^. As a simple further check we have also computed the corrections to 
the masses of bulk particles. In particular we have focused on the n = modes of 
fields with spin s = {0, 1} and bulk mass {mQ,m']^}. Again, we obtain 



where the gauge dependent part is the same as for brane modes but the physically 
meaningful piece A'^ is, as expected, different. Notice that bulk particles do not 
couple directly to branons, so that there is no analogue of diagrams (d) and (e) for 
them. On the other hand the bulk modes couple directly to the hi^ and hij pieces 
of the graviton field, which was not the case for the brane modes. In view of these 
differences, the fact that the gauge depended piece is always the same is a rather 
non trivial check. 

A concluding remark on the terms found in [^, ^ ^, |^ is in order. These 

terms come only from diagram (a), so that the branons play no role in the cancellation 
of these effects in physical quantities. Furthermore, it is clear that terms of this type 
could not be physical, as they cannot arise in the ^ = 1 gauge. However in gauge 
dependent quantities they can appear. In the appendix we give the expression of 
the scalar self-energy at the off-shell point Pext = 0, where these unphysical effects 
are indeed present. Notice that if they appeared in physical quantities there would 
really be an enhancement of the result by some power of the radius R (IR divergences 
cannot be thrown away!). 

So far we have only considered observables that do not depend at tree level on 
the size R of the compact extra dimension. A gauge invariant result is obtained in 
a slightly more complicated way when one considers observables like mo/Mpi or the 

^The gauge dependence of pole masses in quantum gravity was already found and discussed in 
ref. [p7[. In the next section we will give a simple geometrical explanation. 



dm' 



s 



2m'; G + m[ 




(3.5) 
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ratio between pole masses of different KK excitations. The reason is that R itself is 
gauge dependent: a discussion of this issue, including a geometrical explanation of 
this statement, is presented in subsection [O . 



3.2 Gauge independence of physics and geometry 

In this section we want to extend our discussion of gauge invariance to generic ob- 
servables that depend on the size R of the extra dimensions. To be concrete we 
will compare two gauge choices, unitary (U) and de Bonder {DD). The compact 
manifold is assumed to be a ^-torus. 

The previous results could be restated as follows: in order to get the same physics 
in the U and DD gauges, all the tree level mass parameters mu and m£)£) in the two 
gauges should be related by 



[1 + - Gdd)] = ml^Xi [m^] = 2, (3.6) 



where the G's are the universal quantities given in appendix C. Similar relations 
hold for parameters with different mass dimensions. This is equivalent to taking 
as background metrics Xitjmn and timn in respectively the U and DD gauge, but 
keeping the same tree level mass parameters {i.e. moD)- This is easily seen because 
the tree level Klein-Gordon operator in U gauge is Xi^rj^^^ dudN + i^'dd- 

This argument is basically correct, but not completely. The point is that since 
the space is not isotropic (there are 5 compact directions) the metric rescaling factor A 
does not have to be the same for all directions. Then, compatibly with the symmetries 
of the system, we expect in general the backgrounds 

with Ai 7^ A2. These relative backgrounds have to be chosen in order to get the same 
results in the two gauges. Notice that we keep the same periodicity ~ + 2t[R on 
the torus. Then, at tree level, the proper length of the period of the torus is rescaled 
by a factor in the unitary gauge. In the same gauge the mass shell condition 
for the mode {rij} is 

[V^'^d^d, + + Aim2)(/)„ = (3.8) 



SO that R as defined through KK masses is rescaled by \/X2/Xi at tree level, and not 
by y/X2. Finally the tree level ci-dimensional Newton constant Gn = 1/{M^^~'^R^) 
in the U gauge is given by = GnX2^^'^X\ '^^'^ . By writing Aj = 1 + q, at lowest 
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Figure 2: The corrections to the Planck mass is obtained by combining corrections to the 
graviton propagator (we show the relative Feynman diagrams. Diagram b contains a ghost 
loop) with corrections to the graviton/matter vertex and with corrections to the matter 
propagator. 



order in the q we then have the tree level relations 

mu =mDD{l + 

i?c/=i?DD(l + |-|) (3.9) 

{Gr,)u={GN)DD{l + {l-^)ci-^-C2) 

where the radius is here defined through the KK masses. In our calculations so far we 
only considered the masses of brane modes or bulk zero modes, which do not depend 
on the radius at tree level. This is why one universal rescaling Ai was enough to elimi- 
nate spurious gauge effects. By considering the quantum corrections to KK masses or 
to the d-dimensional Newton constant one finds extra gauge dependence. However 
we have checked by explicit calculations that it can all be eliminated consistently 
with eq. ( p.9|) . The corrections to KK masses represent just a direct generalization 
of the computation of the previous section. On the other hand, the Newton constant 
requires to compute also the correction to the graviton-matter vertex and to the 
graviton propagator. Few relevant Feynman diagrams are shown in fig. ^ This is a 



lengthy computation^ upon which the gauge dependence in eq. (3^) is a non-trivial 
check. At one loop, we find 

" (3.10) 

" n 

It is interesting to re-derive the quantities Ai,A2 in a purely geometrical way. For 
instance A2 is fixed by the coordinate independent proper period of the torus. We 
can easily show this for the case 5 = 1, li = 4 (the latter choice being made just to 
simplify the notation). In order to do so we must (arbitrarily) pick a path around the 
compact dimension, and make sure that working in different gauges the path is kept 
unchanged. It is convenient to simply pick the path V defined in unitary coordinates 

^For example in the unitary gauge the diagram |^a is obtained by summing 2.588.740 terms. Ah 
computations in this work have been done with Mathematica ps| . 
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by = 0,y^ = T with r going from to 2ttR. Actually any path in the family 
= const, = r (r = [0,27ri?]) would give the same result as it is equivalent by 
translation invariance of the background''. Notice also that in a general coordinate 
choice, x'^ is not constant along V. In an arbitrary coordinate system the definition 
X*^(r) of V depends also on the metric itself. It is easy to work out this dependence. 
The quantity 



j^\/gMNX''X^dT^ (3.11) 



must be gauge independent, being the expectation value of a gauge invariant oper- 
ator. Comparing the calculation of L in the U and DD gauges, we get at 1-loop 
order 



V {^5fM^5u )dD + -:[ — ^ \%5 ^55 )dD 



4 V ml 



In the unitary gauge only the scalar zero mode (radion) h^^ contributes to L, while 
in the DD gauge extra contributions from KK graviphotons and graviscalars show 
up (the latter is zero in dimensional regularization) . However in the U gauge there is 
the tree level term C2- Notice that in both gauges the gravitational field h is defined 
to have no tadpoles. This equation fixes C2 and the result agrees with what found 
for the physical parameters. 

Finally we can fix ci by considering the volume element in the two gauges 



/ ^ d'y d'^x) =11 y^d'y d^x 



(3.12) 



DD 



where we have fully integrated on the torus . This equality ensures that the non- 
compact coordinates x represent the same physical distance in the two gauges. The 
mass of a particle, as defined by the x dependence of the propagator, has then to be 
the same in the two gauges. Taking the background into account, eq. ( |3.12| ) reads at 
1-loop 



(27ri?)^ 



{2t^RY l + {d)DD 



(3.13) 



^This is an important point since by construction the unitary coordinate frame, defined by the 
request that 555 and be independent of y^, is truly a family of gauges. This is because the 



unitary form of the metric is preserved by the "zero mode" coordinate changes x'^ 



y5 _j, y5 ^_ (^(^x) (corresponding to 4-dimensional diffeomorfisms and to the circle isometry). Then 
since our path V is defined in a family of gauges it truly designates a family of paths. It is manifest 
that this family corresponds to the paths related to V by translation in x. They all have the same 
length. 
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Figure 3: A diagram of order t^. 



where 



= ^ _ 2,^^^^^)/.£U5,-"). (3.14) 



Using the previous result for C2 we determine Ci in agreement with eq. ( |3.1CI| ). 
3.3 Higher orders in the brane tension 

We have only studied the terms of zeroth order in the tension r, but things should 
work out in a similar way order by order in r. These higher order effects come not 
only from branon insertion in the diagrams of fig. (la) and (lb) but also from extra 
tadpoles. Indeed at order r there is already at the tree level the tadpole of fig. (|^). 
It corresponds to the brane self gravitational field. Of course if we treat the brane as 
a thin object this field is infinite at the brane itself. This is a UV divergence which 
is eliminated by adding the suitable counterterms (amounting to a renormalization 
of the unit length on the brane). Applying our regulator (see appendix |^) we get 
from fig. (^, for d = 4 

where Ii is a constant defined in appendix 0. The R dependence is insensitive 
to the UV cut-off: it measures the deformation of the brane self field due to the 
finite volume, so it is a well defined quantity in the EFT approach. Notice also 
that fig. (^) is a brane-to-brane exchange of a bulk graviton like those considered in 
various phenomenological studies 0, At one loop fig. (H) is dressed into extra 



tadpole diagrams: we expect that inclusion of these tadpoles will be essential to get 
gauge independent results at linear order in r. 

3.4 Brane in an orbifold 



In section |2]^ we explained that, for a brane living on a smooth space, the branons 
have to be kept in order to preserve general covariance. Then we have explicitly 
shown that the branons are needed to restore reparametrization gauge independence 
of quantum gravity corrections. In this section we show an example of how things 
work for a brane stuck at a fixed point of an orbifold. Now the brane cannot move, 
i.e. there is no branon degree of freedom. But at the same time the group of 
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difFeomorfisms is also changed. Indeed when deahng with fixed points it is even 
superfluous to talk about a brane: at these points we can localize degrees of freedom 
and interactions respecting the orbifold reparametrization invariance. For simplicity 
we will consider the simplest case of a brane in W'' x S^/Z2. The space 3^/2,2 is a 
line segment, obtained identifying points in a circle of radius R according to the Z2 
reflection: 

y-^27rR-y] y e [0, 27rR] (3.16) 

which has 0, tiR as fixed points. The invariance of the line element rfs^ under Z2 
implies that under the orbifold refiection the metric components /i^^, hij and the 
ghost field f]^ are even, while /i^j and rji are odd. A generic field f{x,y) can be 
Fourier decomposed according to its parity: 



/(a;,2/) = J]/(")(a;)v^„(l/), ^„(|/) = 

n=o \Ksm{ny/R) odd 



Gn cos{ny I R) even 

(3.17) 



\'2-kR y-nR 

Odd fields do not have a zero mode. We can use the same gauge fixing for reparametriza- 
tion invariance as before. 

The group of diffeomorflsms on the orbifold is deflned by the transformations 
x'^ — > f'^{x, y),y ^ y) with even and p odd under orbifold reflection (both 
and p have period 2t^R) . Notice that the boundaries y — Q and y — ttR are left 
fixed. A brane aX y = remains a brane at 7/ = in all reference frames. Even if we 
do not let the brane fiuctuate we still obtain consistent results. On the other hand 
for a brane at a generic y ^ 0, nR, its position depends on the reference frame and 
we are forced to let it fiuctuate. What is special about the fixed points is that we 
have thrown away enough gauge degrees of freedom {g^^{y = 0,nR) = 0) that we 
can live without branons. Let us see this explicitly. 

The computation of the gravitational corrections in the orbifold geometry is 
similar to the previous ones, but with some important differences. Consider the case 
of a brane sitting at a generic point y. Contrary to the circle case, the y dependence 
in the coupling between matter on the brane and gravity does not cancel out in 
physical amplitudes; for instance, the cross section for the production of an individual 
KK graviton mode is proportional to cos^ (ny/R). It is not a surprise that this factor 
depends on y: is not an homogeneous space. Similarly, the branon contribution 

in graviton loops gets multiplied by a factor s\v?{ny/R), showing that their presence 
is not necessary when y = 0. However, having altered by a y-dependent factor the 
relative weight between graviton and branon effects, we apparently no longer get a 
gauge invariant result. We now show that we must take into account a new type of 
graviton tadpoles that were absent on a homogenous space. 
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The conservation law of fifth dimensional momentum is altered since some of the 
harmonics are projected out by the Z2 symmetry. In a vertex with three lines carrying 
momenta Ui > 0, i = 1,2,3 along the fifth dimension it reads ni ± 77,2 ± ns = 0. The 
propagator of matter on the brane is corrected by new tadpole diagrams, like (Ic), 
but with non zero extra-dimensional momentum 2n on the tadpole graviton line. The 
blob in fig. (Ic) can be either a graviton loop or a gravitational ghost loop. Notice 
that, while tadpoles with a zero momentum internal line (n = 0) are assumed to be 
exactly canceled by a suitable stabilization mechanism, the same type of diagrams 
with non zero n must be taken into account and they are crucial to recover gauge 
invariance for brane observables. Let us focus for instance on the mass correction 
Sitlq for a scalar on a brane at a generic y. Notice that for y ^ 0, ttR the brane is 
free to move, so branons must be kept. We get 

Sml = ^ [2 cos^ (ny/R) ^("^(gravitons) + 2 sin^ (ny/R) F^^^ (branons) + 

n 

cos{2ny/R) F^") (tadpole)] = ml (2G{MdR) + Ao(mo, y, R, Md, fi)^ 

(3.18) 

The contribution from the nth KK mode in the graviton diagrams of fig. (la,b) is 
exactly the same as in the torus, except for an overall factor 2cos'^{ny/R) coming 
from the graviton wave function. The contribution F*^") (branons) from diagrams in 
fig. (ld,e) gets instead an overall y-dependent factor 2 sin"^ {ny / R) . The final result 
for 6ml has the same structure of eq. (p.4|), but now the gauge invariant piece Aq is 
a function of y. Finally, the tadpole contribution comes in the right way to cancel 
mode by mode the ^/-dependence in the gauge variant term G. This is consistent with 
the mass correction for the zero mode of a scalar propagating in the bulk, which has 
the same form as in the torus case (see eq. ( |3.5D ), and it represents a non trivial 
check on the result. In the special case of a brane sitting at the fixed points 0, vri? 
the branon contribution vanishes and gauge invariance of the pole mass is met just 
through tadpole diagrams. 

As a final remark, we notice that because of the modified momentum conserva- 
tion law in the orbifold, the zero mode of a bulk field mixes with its KK excitations 
at one loop level. The relevant diagrams are those in fig.s (la,b,c) with discrete 
momentum n in the internal loop and 0, 2n in the external legs. This effect however 
is relevant only at order k'^ and can be safely neglected. 

4. Phenomenology 

In the previous sections we have explained how to consistently compute quantum 
gravity corrections using an effective field theory (EFT). A possible physical ap- 
plication is the computation of graviton loop corrections to electroweak precision 
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observables (EWPO) and to the anomalous moment of the muon in brane models 
with large extra dimensions and a TeV-scale D-dimensional Planck mass. Unfor- 
tunately our knowledge of the low energy effective theory of gravity only allows to 
reliably compute corrections of little phenomenological interest. Basically, the EFT 
allows to compute those contributions that are saturated in the infrared, i.e. at the 
scale of the relevant external momenta. For instance, the calculable corrections to 
EWPO go like {Mz/MDf+^ (or {Mz/Mof+^lriMz) by simple dimensional analy- 
sis. These effects go to zero very quickly when M/j is raised, becoming negligible 
already for Md below a TeV. On the other hand, the contributions from the region 
of large virtual loop momenta gives in principle a much larger effect. However, being 
saturated in the UV region, where we do not control the EFT, these contributions 
are not calculable. This problem already affects tree level virtual graviton effects. 

We can however estimate graviton effects by introducing an explicit UV cutoff A. 
The corrections to EWPO will scale like M'^A.^ / M"^^ . The unknown physical cutoff 
could perhaps be produced by string theory, or could be related to the inverse brane 
width or even to just the brane tension ^0[. Since we do not know we must keep 
A as a phenomenological parameter and discuss its physical meaning and plausible 
value. 

Virtual graviton corrections (even at tree level) cannot be computed from Ein- 
stein gravity as much as electroweak quantum corrections cannot be computed from 
Fermi theory. In the latter case the complete theory is known and perturbative: by 
comparing to the full theory one sees that correct estimates are obtained by cutting 
off power divergent four-fermion loops at a "small" scale A ~ M\y ~ gG-p rather 

— 1/2 

than at the larger A Gp . At least at a qualitative level, the gravitational A can 
be given a similar physical meaning. 

4.1 Strong vs weak gravity: NDA estimates 

Therefore we first identify the value A^ of A that corresponds to strongly coupled 
quantum gravity®. This can be done by adapting to our case the naive dimensional 
analysis (NDA) technique developed to estimate pion interactions below the QCD 
scale I^TI (NDA has already been applied to brane models [0). NDA allows to 
estimate the size of the effects from a strongly coupled theory up to coefficients of 
order 1 but including all the geometric dependence on powers of vr. By applying 
NDA, we estimate 

A|+^ ^ 7r2-^/2r(2 + 5/2)Ml+^ (4.1) 

In the range of interesting 5, A5 is not much larger than M^. 

We first discuss the particular case A = A5: diagrams with any number of 
graviton lines give comparable contributions, and NDA allows to estimate their size. 

^In the context of string theory this corresponds to a situation where the string coupUng is 
essentially at the self dual point. 
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Tree level exchange of gravitons generates the effective dimension 8 operator T = 
T^^ — T^^/ (5 + 2) 0, ^ 0. Its coefficient in the effective Lagrangian is divergent, and 
NDA estimates it to be ~ 7r^/A|. This operator is however not the most important 
in low energy phenomenology, because at loop level gravitons generate dimension 6 
four fermion operators with coefficient ~ 7r^/A|. On the other hand the operator 
W^j^B^'^ Wt'^H is generated with coefficient ~ (y'2fi'i/A| with no vr^ enhancement. 
This property is shared by other operators that require the exchange of virtual gauge 
bosons. This is because we are assuming that the weak gauge couplings remain weak 
up to the cutoff. 

By drawing a few Feynman graphs one can see that tree level exchange of gravi- 
tons (and therefore the operator T) does not affect precision observables at the 
Z-resonance. Moreover the four fermion operators induced by double graviton ex- 
change are of neutral current type, so they do not directly affect /i decay and are 
therefore not constrained by high precision data, /i-decay is affected by one loop 
diagrams with a W and a graviton: their coefficient is only k, g\j t^^. 

Indeed by a simple analysis one finds that all dimension six operators that affect 
EWPO have a coefficient ^ g\l ~ V^l- As shown in Q EWPO set a bound 
A5 > (5 ^ f 0) TeV on a generic set of dimension 6 operators that conserve baryon, 
lepton and flavor numbers and CP. This bound seems rather strong when compared to 
the sensitivity to direct graviton emission expected at the next colliders. Furthermore 
since m^op ~ 175 GeV a real solution of the hierarchy problem should cutoff the 
quadratically divergent top correction to the Higgs mass at a much lower value of 
A ^ 300 GeV. Our assumption A = A5 corresponds however to one of the most 
constrained scenarios: LEP data strongly disfavor new strongly coupled physics in 
the electroweak sector. The situation becomes worse if we assume that also the gauge 
couplings get strong at A = ks- 
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In order to obtain a more acceptable phenomenology one can assume that the 
UV cut-off A happens to be smaller than As, so that gravity does not become strong 
and dominant graviton corrections to EWPO are dominated by one loop diagrams 
(presumably a complete theory will not contain only gravitons) . In the next subsec- 
tions we 'compute' the graviton corrections to the electroweak observables (expressed 
in terms of the ei, €2, £3 parameters |12[) and to the anomalous magnetic moment of 
the muon. In agreement with NDA estimates, the final result is of the form 

where the factors of order one depend on the choice of cutoff. Not knowing which is 
the physical cutoff, we use dimensional regularization: with this choice loop integrals 
do not give powers of A. However, since we are considering a higher dimensional 
theory, powers of A arise from divergent sums over the KK levels of the gravitons. 
A different choice of the cutoff would give different results. 

In fig. ^ we summarize the present situation of collider graviton phenomenology 
by collecting the various bounds in the plane {Md,A/Md): 

• The vertical bound comes from emission of real gravitons at LEP2 and 
Tevatron |Q . It does not depend on A (as long as the energy of the collider is 
less than A) because it is the only bound on really computable effects. 

• Virtual exchange of gravitons at tree level generates the operator T. Its coeffi- 
cient depends on the cutoff A so that it cannot be computed from the low energy 
EFT (in the literature there exists a variety of estimates 0, ^, ^ , freely dubbed 
"formalisms", and a corresponding variety of experimental bounds [^). The 
coefficient can be estimated to be ^ tt^5{5 + 2)K^~'^/2{5 - 2)A|^^ The exper- 
imental constraints give the slightly oblique bound in fig. |. 



• At one loop gravitons affect precision observables and in a way that again 
depends on the cutoff. The green line shows the values necessary to produce 
the observed anomaly in a^. The bound parallel to it comes from precision 
observables. 

If the cutoff A is due to quantum gravity, K/Md parameterizes how strongly coupled 
gravity is: this explains why virtual graviton effects give the strongest (weakest) 
bound when A>M£) {A<M£)). Strongly coupled gravity is obtained for A/Md ~ 
(1 ^4) if (5 = 3 and for A/Md ~ (1 ^2) if 5 = 6. EWPO bounds have been estimated 
in a conservative way, assuming a typical 0.1% error. We see that setting A = Md 
as assumed in many analyses is a significant but arbitrary restriction: A is a relevant 
free parameter. In the most generic case the cutoff could even be not universal, so 
that different corrections are cut off by different A. We repeat that bounds that 
depend on A can at best be considered as semi- quantitative. 
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The presence of a cut-off A can have an impact on the studies of graviton emission 
at future colliders. If A is smaller than a/s, real graviton signals are suppressed (but 
some new physics should show up). On the other hand, if A/Md is too big, real 
graviton signals (7+ missing energy) are forbidden by precision tests or subdominant 
with respect to 7+ missing energy effects due to dimension six operators like eevv [^], 
generated by virtual gravitons at one loop with coefficient ~ 7r^A^'^+^/A|'^'''^. However 
there exists a range of Mr, and A (not too small and not too large) where real graviton 
emission is the dominant discovery mode. For instance one can see this by considering 
the case of a e+e~ collider at = 1 TeV @]. 



Can the apparent excess a^^P — = (4.3 ± 1.6) ■ 10~^ recently measured by 
be produced by gravitons without conflicting with the EWPO bounds? In the SM, 
electroweak corrections have been clearly seen in the e^, but only affect at a level 
comparable to its present experimental error. The naive (and maybe correct) expec- 
tation is that even in the gravitational case the are a more significant probe than 
a^. However, taking into account that we can only perform estimates, it could not 
be impossible that the anomaly in [^] be produced by gravity without conflict- 
ing with the EWPO bounds, even if the physical cutoff has a 'universal' nature (for 
example if it is related to the size of the brane) as assumed in fig. ^. If this is the 
case, improved measurements of the parameters should be able to find a positive 
signal. 

4.2 Electroweak precision observables 

As discussed in the previous sections, unphysically large corrections cancel out when 
correctly computing physical observables. Previous analyses have studied certain 
combinations of the vacuum polarizations of the vector bosons 

n;f.(^') = -lV^..^'Kk^)+k,K terms, i,j = {W,Z,^} 



known as S", T, \J parameters |]TT| , often employed to parameterize new physics present 
only in the vector boson sector. However these are not physical observables because 
gravity does not couple only to vector bosons^. As found in 0, in the unitary 
gauge gravitons give corrections to such parameters that unphysically increases with 
increasing Mb- 

Since graviton loops are flavour universal (and neglecting the bottom quark mass) 
gravitational corrections to the various EWPO can be condensed in three parameters 
that are usually chosen to be 61,62,63. The corrections to the physical EWPO are 
obtained by combining in a non immediate but standard way []T^ various form factors. 
Specializing the general expressions to the case of gravity, the 6 parameters are given 



^Various studies on different new physics scenarios use the S", T, U approximation outside its 
domain of appUcabihty. 
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by 



^1 =^5g - ^ - ^ - n'^z(M|) (4.2) 

z 

e,=2c=^«-f -^+.^^-c^n^,(M|) (4.3) 

£3 =2^*9 - c'- - c'ni2(M|) (4.4) 

a 



where 



• 5Ml = —Ilii{Mf) are the correction to the pole mass of the vector bosons and 

• Sa — — n^-y(O) is the correction to the electric charge; 

• 5g is the common correction to the vector and axial form factors (gravity re- 
spects parity) in the Z^^ff interactions of an on-shell Z boson 

g _ 

-«^/7m(^v - 755a)(1 + Sg)f 
excluding the contribution from the Z vacuum polarization. 

• 5G is the correction to the /i — > ePg^^/i decay amplitude. 

• = 1 _ c2 = [1 - ^l-A7ra/V^GM^]/2 ^ 0.2311 

Although it would be straightforward to perform a complete analysis, we will only 
study the gravitational correction to the combination 

chosen because it only involves the simplest-to-compute form factors. Physically, 
this observable amounts to testing the tree level SM prediction — cMz using 
the value of the weak angle given by the forward-backward asymmetries in Z ^ 
decays, not affected by graviton loop effects. The experimental value of e (obtained 
from a fit of LEP and SLD data) is e = (12.5 ± 1)10^"^ and agrees with the SM 
prediction (for a light higgs). The gravitational correction is given in appendix C 
in terms of Passarino-Veltman functions. Since the heaviest KK give the dominant 
effect, we can explicitly write the graviton effect in the limit m„ ^ Mz as 



M|(47r) 



40 + 25(5,1 . jl\ 424 + 5465 + 137^2 



(-+1-37) + 



6 + 35 ml^ 18(2 + 5)^ 
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where we set (i = 4 — 2e. We can estimate the graviton correction by keeping only 
the logarithmic term, setting = A and cutting off the sum at n < RA. This gives 

g^M| / A \ ^ 5(8 + 5^) 
Ml \Md) 48r(2 + 5/2)vr2-5/2- 

The result has a strong dependence on A and the numerical coefficient is specific of 
the form of the cutoff that we have chosen to employ. We explicitly see that spurious 
IR divergences do not affect this physical observable. 

Notice that e is suppressed by a power of but it is not affected by theoretical 
uncertainties in 6a. Therefore precision searches for Md could be improved by a 
factor ~ 3 if, by producing ~ 10^ Z bosons at an ee linear collider, the errors on 
Mw and on the effective weak angle extracted from the leptonic asymmetries could 
be reduced by a factor ~ 10. 

4.3 Anomalous magnetic moment of the /i 

Since the fi anomalous magnetic moment is zero at tree level, the reparametrization 
gauge dependence of the unit of mass does not affect the one loop gravitational 



correction to a^. Only few Feynman diagrams contribute. As noticed in the 1/e 
poles cancel out when computing the loop integrals using dimensional regularization 
around d = 4. At leading order in we find, again using a sharp cut off for the 
sum over KK modes 



Sa,, 



ml f A\^ 34 + 11(5 



MI\MdJ 961(2 + (5/2)7r2-V2- 

Apparently this result agrees with the one found by |]Tn|^°- Again the result strongly 
depends on the value of the cutoff A. We cannot claim that it has the same sign as 
the apparent excess recently measured by we have employed dimensional reg- 



ularization for loop integrals but other regularizations (e.g. dimensional reduction, 
Pauli-Villars, ... ) would give a different result. Unlike e, 5a ^ is a sum of contri- 
butions from graphs with different graviton interactions. One can obtain any sign 
for e.g. by cutting off /i/i/i and 77/1 vertices with different form factors: 5a^ is 
finite but dominated by loop momenta around the cutoff. In particular one gets, at 
leading order in 

Sa^ = 

if the cutoff acts in the same way on both type of contributions. This is for example 
the case of a Pauli-Villars cutoff on the graviton. By working in the De Bonder 
gauge (where the only dependence on the graviton mass comes from the l/(fc^ 
factor in the graviton propagator) and knowing that is dimensionless and finite, 
it is not difficult to realize that it is zero. 



^"pOt separately computes the gauge-dependent 'graviton' and 'radion' contributions in the uni- 
tary gauge. We find a different result in both cases (the radion coupling used in is valid only 
on shell), but this discrepancy luckily cancels out when summing the two contributions. 
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5. Conclusions 



We discussed various subtleties that arise when computing quantum gravity 1-loop 
effects in models with large extra dimensions and matter confined to a brane. Our 
computations are based on an effective field theory (EFT) description of quantum 
gravity and of the brane. A sensible result is obtained after correctly identifying phys- 
ical observables and after taking brane fluctuations into account. Graviton tadpoles 
are relevant for branes living in non-homogeneous spaces (like orbifolds). For branes 
living at orbifold fixed points consistency is met, as expected, even in the absence 
of brane fluctuations. In particular we explain in a geometric way why the units of 
length in 'longitudinal' and 'transverse' directions depend on the reparametrization 
gauge fixing procedure. 

We regard these results as theoretically interesting, although the truly calcula- 
ble effects in the EFT approach have a limited phenomenological relevance. The 
most relevant effects come from the region of large virtual momenta where the EFT 
description breaks down. This is why in the second part of the paper we have 
abandoned the strict EFT approach and modeled these UV effects by introducing 
a hard momentum cut-off A. This is the best that can be done, without having 
a fundamental theory that allows real computations. We stress however that our 
previous understanding of how to get gauge independent results is still important in 
this phenomenological approach. As an application we have studied virtual graviton 
corrections to precision observables and to the muon anomalous magnetic moment, 
focusing on models with large extra dimensions. Even at tree level, virtual graviton 
effects are divergent and must be regulated. Virtual graviton effects in collider phe- 
nomenology have been so far studied assuming a particular value of A. However A is 
an important free parameter that — at least at an qualitative level — controls how 
strongly coupled gravity is. Depending on the value of A, one loop effects can give 
the dominant bound on low scale quantum gravity. 

Acknowledgments 

We would like to thank Ignatios Antoniadis, Gian Giudice and Angel Uranga for 
interesting discussions. R.C. is indebted with G.Policastro, R.Sturani and especially 
with A.Gambassi, M.Gubinelli and A.Montanari for many useful discussions. This 
work is partially supported by the EC under TMR contract HPRN-CT-2000-00148. 

A. Graviton propagators in the unitary gauge 

We derive here the propagators for the physical fields in the unitary gauge. Expand- 
ing the metric qmn = Vmn + i^^mn around the flat space solution we obtain the 
Lagrangian in eq. ( |2.1| ). We recall that we decompose the D dimensional graviton 
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as h^y, h^i and hij (where /x, v are 4-dimensional indices and z, j span the extra 
5 dimensions). In this appendix we fix for simphcity d = A. Due to D-dimensional 
reparametrization invariance not all the components of these fields correspond to 
propagating degrees of freedom. The physical fields are the ones contained in the 
Riemann tensor 2Rrmsn = h(^FiN,MS) — h(MN,RS) (as in the electromagnetic case the 
physical fields are contained in the field strength tensor) 

Gfj,u = — '^didjRi^ju = hfj_y — did(^^hy)i + d^dydidjhij (A.l) 
Vfj_i = + 2djdnRjf_iin = h^i — didnh^n — d^djhij + d^didjdnhjn (^-2) 

Sij 2,0jYi0nRimjn ^ij 9„0(^ihj'^^^ -\- OiOjO^Onhfyifi (A. 3) 

For simplicity, the above equations are written assuming units such that didi = 1. 
These expressions can be written in a compact form by defining = dih^i, P = 
didjhij and Pi = djhij — diP, where di = dij ^fdjdj. These considerations suggest to 
rewrite the Lagrangian in terms of a new set of fields G^,^, V^j, S'jj, Q^, Pi, P 
related to Hmn by 

V =G>.u - + -y^)H + d^Q, + d,Q^ - d,d,P (A.4) 

h^j =S,j + ^(%- - ^^^J)H + d^Pj + d,P, + d,d,P (A.5) 

Ki =V^i + d^Pi + (A.6) 
and subject to the constraints 

^^V^, = d,S,, = d,P, = 0, S,, = 

We have introduced the field H in order to make Sij traceless. By choosing = 
3(5—1) / (5+2), H is canonically normalized. The D-dimensional Lagrangian becomes 

£ = - lH{n + dl)H - \s^^{p + dl)S,, - V^%u + dl)^,. - d^d^^V^ + 

2 2 (A.7) 

-G%u + dl)Gl - -G^^(n + dl)G^, + G^^d^d^G^ - G^^d^d^G^r 

As expected, it does not depend on Q^, Pi and P and there is no mixing between 
the fields G^y, H , Sij, V^,. It is now trivial to perform a mode expansion: the extra 
dimensional Laplacian dl becomes a mass term. The propagators can be obtained 
by inverting the kinetic terms in eq. ( |A.7| ). It is useful to show explicitly how the 
'graviton' G^^ and the 'scalar' H combine to give a unitary-gauge propagator equal 
to the de-Donder propagator, up to longitudinal terms. For example, the h^^h^/y' 
propagator is 



I'^'n r) ' /' 2 



I,, I 
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In particular we see that we cannot omit the 'scalar' contributions, if we want to 
obtain a gauge invariant result. It would be easy to include a small mass term for 
the H^"^^ fields, eventually generated by the unknown mechanism that stabilizes the 
size of the extra dimensions. 



B. Graviton vertices 



We define (7^1. = "q^^ + nh^^, g = \detg^^\ and give explicit expressions for the 
expansion up to second order in the graviton field h^^, of 



^g'^'^ =1 + KB^'^^^h^p + K'B'^''''''^'Kph,s + • • • 
In the vierbein formalism the spin connection is given by 



] '/ma 5 



(B.l) 
(B.2) 
(B.3) 



(B.4) 



We expand around the flat background 5^, = 5^+ k 6^. As discussed in section 
the gauge choice = allows to express b^^ in terms of /i^i. 



fJ,U 



(B.5) 



Using 



3!) and 



(B.6) 



one can find the gravitational couplings for fermions. 

With these expressions it is straightforward to find the graviton vertices arising 
from Lagrangians like 



9 



-oPPo'^^F F 



m 
~2 



^g^'A.A, + 



where F^, = d^A, - d,A^ and D^^ = d^^ + icu/^7,bV, 7^' = i[7",7l- The 
expansion in powers of h is easily obtained from 



^1 Tr/i 2 
=1 + + « 



Tr^ h Tr h"^ 



Tr^ h Tr /i Tr Tr /i^ 
^ + 



48 



6 



+ 



(B.7) 
(B.8) 
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Therefore 



a/3 



^pVpCFCxP 

^ipvpaafS-yS 



(B.9) 

(B.IO) 
(B.ll) 
(B.12) 



L4 



(B.13) 



These expressions are vahd in any number of dimensions. Brane fluctuations can be 
incorporated in h^^, as discussed in eq. ( p.l2|) . 

To compute the corrections to the graviton propagator and to the graviton vertex 
it is necessary to have the 3 and 4 graviton interactions. They can be easily derived 
by expanding the Einstein-Hilbert Lagrangian in powers of the graviton field using 



e.g. Mathematica [2S]. For this reason we do not write explicitlt the long expressions 
for such vertices. 



C. Results 



In this appendix we collect the explicit results for the corrections to the propagator 
of a spin 0,1 particle confined on a brane with dimension d living in R'^ X As 
discussed in sec. 3, generically the correction to a physical quantity O with canonical 
dimension do has the form 



bOjO = do G{MdR) + A(0, R, Md, ^J) 



(C.l) 



where the gauge dependence is encoded in the function G. The splitting in a 'gauge- 
dependent' and 'gauge-independent' part is ambiguous unless a reference gauge is 
chosen in which by definition one sets Gref = 0. We choose in the de Bonder gauge 
Gde Dondcr = 0. All the rcsults for physical quantities are computed in this gauge. 
The results are expressed in terms of Passarino-Veltman functions Aq, -Bq,!, defined 
in appendix 0. 

For the pole mass correction for a scalar (s = 0) and a massive vector (s = 1) 
on the brane we find 



G 



unitary 



V4 ,. Ad' + {5-^)5 + d{25-l)\ 



(C.2) 
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8d[mlid - 2) - 2m2(d + 5-3)] Boiml ml, ml) + 2d5{d - 2)mlBi{ml, ml, ml 

(C.3) 



Ai = ^-^^ > <^ f2(d, 6)Ao(mi) + (d - l)Md, 6)Ao(m' 

647Tm^-^d{d-l){d + 6-2)^^^ ^ ' ; ov i; v jjn , J ok n 



n 



8d[2ml{d - + 5 - 3) + ml{d -2){d + 26 + l)]Bo{ml, ml, mj) + 
2d{d - 2) [2d^ + 3d{5 - 2) - 75] m^5i {mj, ml, mj) 



(C.4) 



with 

f\{d, S)=4 [2d^ - d{6 - 6) + 2d\6 - 3) + 26{6 - 2)] 
/2(rf, 6) = 2d{2 -d)[A + 2d^ + 5 + d{35 - 2 



(C.5) 

We computed also the graviton correction to the photon propagator, verifying that 
it is transverse if one uses the simple gauge fixing of eq. ( p.l9|) . If instead the gauge 
fixing function contains the graviton field, in general transversality will be lost, due 
to a modification of the related Ward identity (of course this does not mean that the 
photon acquires a mass). As the simple QED case. Ward identities imply that the 
photon self energy at zero momentum gives the correction to the electric charge. We 
find 

^ , , ,^ {d-4)[d^ + d\S-5) + d{8 - 36) + 26{6 + 2)] 

Ag = — > An m„ ^-r-^ (C.6) 

^ ^ 327rm^-^d{d + 6-2) ^ ' 

For comparison with the existing literature, we also write the expression of the scalar 
self-energy at zero momentum S(0). As discussed in the text, for this unphysical 
quantity the gauge dependent part is non-universal. In the de Bonder and in the 
unitary gauge we find 

mi 



S(0)de Donder „ Md-2u , A _ 5Z 



327im^-\d + (5 - 2) ^ K - ml) ^^^^^ 
[mlg^{d, 6) + mlg2{d, 6)] Ao{ml) + 2d{6 - 2)mlAo{ml) 



= ^ X V 



327r2M|-"(d + 5-2) 



„ m4(m2-m2) 



2ml [Kd{^ - 2) - 2mlml{6 - 2) + m^(d + 6-3)] AqK) + (C- 
{d-l)m1[ml{d'' + d{5-2)-2) -ml{d-2){d + 6)]Aoi 
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with 

gi{d, 5) = d\d -1)-Ad + S{A + d){d-l) + 25^ 

g^(d, 6)=46-{d + 6) [26 + d{d - 1)] ^^'^^ 

From these expressions it is clear that l/rn^ terms found in [0, |], ^ are an artifact 
of the unitary gauge and have no physical meaning. 

D. Regularized sums and integrals 

The results of our 1-loop computations can be expressed as the sum over KK modes of 
basic Passarino-Veltman functions Generically these expressions are divergent 
and need to be regulated. After doing that one can extract the calculable finite parts 
that are determined by the EFT |jl3|. These are terms that either depend on the 
radius R or depend non-analytically on the kinematic variables. In this appendix we 
focus for illustration on these calculable terms and disregard the uncalculable UV 
saturated contribution, which were the subject of our phenomenological discussion. 

The main point is to regularize the integral and the series consistently; we choose 
for this the dimensional technique, extending the physical dimension of the extra 
space and of the brane 5, d, to generic values 

5 = 5 — e] d = d — e 

rescaling the Planck mass in the Lagrangian as M'j^^"'^ M^'^''^ and 
consequently 

1 _ 1 fi^' 



and taking the limit e — at the end. Defining the Passarino-Veltman functions Aq, 
-Bo,i 



d'^q 



{2tiY g2 _ m 



2 



^ / 9 , ^9 9x —iiAnY f d'^q p ■ q 



p2 J (27r)^(g2 -M2)[(g + p)2 -m2] 

the following expressions are sufficient to compute the gravitational corrections in 
appendix |^ 

Y^M^l) Y.'^n"Bo,iim^,ml,m^) a = 0,1 (D.2) 
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together with the series 

To illustrate the technique we compute explicitly J2n ^oi^"^ ^ ^'^) ■ First of all 
we introduce a Feynman parameter x, rescale the integration variable q q/R and 
isolate the zero point in the series 



neZ^ 



{2nY-^ Jo ^ J ^[g2-n2-a2(x)]2 +y [q^ - {1 - xym^ 

neZ'^ — {0} 

(D.4) 

where a^(x) = i?^m^(l — Then we Wick-rotate and evaluate the first term 

using the Schwinger's proper time method 



y f d^q- / = y r dt ! d'q-^e-'^^'^-'^-'^ 

POO 

=7^2 / ^„r«^"^A _ ll.-!/^«%,l-'^72 



;"00 
^0 



(D.5) 



where we have performed the gaussian integral over q and introduced the special 
function 



oo 



B{s) = J2 e"""'' 

n=— oo 

The integral in eq. (|D.5| ) converges at y ^ oo thanks to the exponential behavior of 
the B function, but it diverges ai y ^ 0. To extract the singularity it is useful the 
property 

B{s) = s~^/^B{-) (D.7) 
s 

which is easily derived from the Poisson formula. Using eq. ( p.7| ) we can split the 
integration interval and change variable y ^ 1/y in the first integral 

/oo _ 
dy [B\y) - 1] [y^-d/^e^y^'"' + yi^^'^-^)!'^^-^-^ N) + (D.8) 

' dy e-™VV/2-3^^572_i) 
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Again, the first integral is convergent, while the second term must be (dimensionally) 
regularized. Because a{x) > in x G [0, 1] and noting that for /3 > 



dy y'^e-l^ly = /3"+^r(-a - 1) 



dy y^+'e 



a+2-l3y 



(D.9) 



we can isolate the divergent piece in the F function through an analytical continuation 
in the physical region d + 6 > 



I 



dy [B\y) - 1] (yi-'^72g-?/™(x)2 ^ ^(j+5-)/2-3g-™(x)Vy) + 



:i-x) 



X 



{d+<5)/2-2 



:i-x) 



X 



d/2-2 



d, 

r(2--) 



(D.IO) 



It's not difficult to verify that the last (divergent) term is exactly canceled by the 
the zero mode contribution of the series in eq. ( p.4|) . Putting together the remaining 
terms we obtain 

-1 



1 



{27tR) 



d-i 



dx 



fi{mR, X, d, 6) + 

where we have defined 
fi{mR, X, d, 6) = 



TTR^m^ 



:i-x) 



X 



2-1 {d+5)/2-2 



r(2-^)} 



(D.ll) 



l-d/2^-yTT0?(x) _^ ,,(c^+5)/2-3^-7ra2(x)/y^| 
yna''(x)yl-d/2^^ _ j/-'5/2y 

(D.12) 

The r function in eq. ( p. Ill) has poles for negative integer arguments and before 
taking the limit e — >■ we must distinguish the two cases of even and odd {d + 5). If 
{d + 5) is even, a logarithmic term appears 



J^d+5~4 



fi{mR, X, d, 6) + 



m 



n(l — x] 



2-1 (d+5)/2-2 



- + log ^ + log ^ - 7i? + 



T{{d+6)/2-r 

1 dF 
F{d, 5) rf7 



X 



e=0 



(D.13) 



34 



where F{d, 6) is a generic function of d, 6 which multiphes the integral in the physical 
amplitudes and the factor comes from the graviton wave function normalization. 
By subtracting just the pole 1/e we get the loop correction in the MS scheme. Notice 
that the finite part contains a scheme independent log m term. When {d + 6) is odd 
we find instead a finite result 



1 1 i\ 

m , m^, m ) 



F{d,6) 



d-4, 



dxx^/''-'' 



^ /i (mi?, X, d, 5) + m'^+^-^ r f 2 - 



7r(l - xf 



X 



{d+5)/2-2 



(D.14) 



Although there is no logarithm, the term represents a scheme independent 

finite effects as it depends non analytically on the Lagrangian parameter m?. The 
same technique can be used to compute the finite part of the other integrals in 
eq. ( p.2| ) and the series in eq. ( p.3| ); here we collect only the final results omitting 
the derivation (for d > 2, a < 5/2) 



-A-K 



{27rR) 



d-2 



dy [B\y)-l]y'/'{l+y'^'~') + 



d-2 d+5-2 
(D.15) 



T 



IT 



Via] 



[ dy[B'{y)-l]{y'^-' + y 



a-l I (5/2-l-o\ 



1 



1 



a S/2 — a 



(D.16) 



Notice that the sum of the Ao(m^) function has no 1/e pole. The special cases Xq, 
are needed respectively to evaluate terms J2 ^o("^o i) results of the previous 

section and the series of eq. (|3.15|) 



Xo = -l 



/oo r 
dy [B\y) - m + y'/'-') - j^] 



(D.17) 



Finally, for [d + S) even 



2 2 2^ 



1 TTl 

f2{mR,x,d,6) + 



d+5-2 



Jld+5~2 



T{{d + 6)/2) 



7l(l — x) 



X 



{d+5)/2-l 



1 , /i^ , 2v/i 1 dF 



1 

-o~6 



(D.18) 
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while for {d + 6) odd 



d + 6 



nil — x) 



(d+<5)/2-l 



(D.19) 



where i = 0, 1 and mo(x) = 1, u\{x) = (a; — 1) and we have defined 

/■oo 



(D.20) 
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